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PREFACE 


The High-Speed Research Program sponsored the NASA High-Speed Research Program 
Aerodynamic Performance Review on February 8-12, 1999 in Anaheim, California. The review 
was designed to bring together NASA and industry High-Speed Civil Transport (HSCT) 
Aerodynamic Performance technology development participants in areas of: Configuration 
Aerodynamics (transonic and supersonic cruise drag prediction and minimization) and High-Lift. 
The review objectives were to: (1) report the progress and status of HSCT aerodynamic 
performance technology development; (2) disseminate this technology within the appropriate 
technical communities; and (3) promote synergy among the scientist and engineers working HSCT 
aerodynamics. The HSR AP Technical Review was held simultaneously with the annual review of 
the following airframe technology areas: Materials and Structures, Environmental Impact, Flight 
Deck, and Technology Integration. Thus, a fourth objective of the Review was to promote synergy 
between the Aerodynamic Performance technology area and the other technology areas within the 
airframe element of the HSR Program. 


The work performed in the Configuration Aerodynamics (CA) element of the High-Speed 
Research Program during 1998 was presented in the following sessions: 


Propulsion Integration 

Analysis Methods 

Design Optimization 

Testing 
The work performed in the High Lift (HL) element of the High-Speed Research Program during 
1998 was presented in the following sessions: 


High-Lift Configuration Development 
Tools and Methods Development 


The proceedings for the Aerodynamic Performance Annual Review are published in three 
volumes: 


Volume I, Parts 1 and 2 Configuration Aerodynamics 
Volume II, Parts 1 and 2 High Lift 


AP Review Chairperson: David Hahne 
NASA Langley Research Center 
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Configuration Aerodynamics Technoloc 
Session 4: Testing 


Goals Demonstrate Significant L/Dmax Gains 


Objectives Robust Analysis / Realistic Aerodynamic Efficient Engine / 
Testing Methods Design Optimization Airframe integration 


Challenges [Watdaton] [Viscous Etec] Multi-Point Conditions | | Power Effects 


Aaniaaehes Analytic Methods Design Fest Programs 
PP and Applications Development and Techniques 


Methods Down Select Nacelle / Diverter WT Database 
Design Integration 


Program Viscous Drag Prediction WT Data Corrections 
——— Tech. Baseline 
Cruise Point Optimization Development High Re. No. Testing 


Multi-Point Optimization Aero S&C PIE Test Program 
Development 


S&C CFD Predictions 





This work was done to support the the WT Database and WT Data 
Corrections programs under the Test Program and Techniques approach for the 
Configuration Aerodynamics Technology Development. 
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Background 
High Speed Research Program seeks to achieve 
new levels of precision and accuracy 
in supersonic performance testing 


AMGoodsell/RAKennelly Slide 3 





Background: 


One of the goals of the High Speed Research Program is to achieve new 
levels of precision and accuracy, particularly in supersonic performance 
testing. 


The picture shown below is that of the 1.675% Baseline TCA Model 2b. 
This is our baseline performance model which has been tested numerous times 
in the past couple of years. 
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» Repeat measurements 
of C, & Cp to reduce 
statistical uncertainty 


. . 
rf * 
. 1 


“GrUISSTi COSTICIENT ~ 


Least-squares 

: polynomial fit used to 
: represent data 
Primarily interested in 
drag at cruise lift 
coefficient 


0.2 
000 0.01 0.02 
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Background (continued): 


As part of the goal to obtain increased precision and accuracy, it 1s 
important to reduce the statistical uncertainty in the wind tunnel data. This is 
achieved by taking between 3-6 repeat measurements of lift and drag over a 
range of angles-of-attack. 


The plot on the left shows an example of a drag polar for the 1.675% 
Baseline TCA Model. The squares are the actual lift and drag measurements. 
The data are represented by fitting a least-squares polynomial curve through 
the points. 


However, we are primarily interested in determining the drag at cruise lift 
conditions. 
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Motivation 


¢ Why do we care so much about drag? 


L drag count (0.0001) ~ 6,000 Ib. take-off weight 


¢ Ultimate program goal is to measure drag 
to within +0.5 count with 80% confidence 
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Motivation: 
Why do we care so much about this cruise drag? 


The reason is that, at supersonic cruise, each drag count costs 
approximately 6000 lb in take-off weight, which is equivalent to about 30 
passengers. Therefore, there is a significant economic incentive to both know 
what the drag is and to try to reduce it. 


The ultimate goal of the Configuration Aerodynamics Test Team is to 
obtain drag to within +/- 0.5 counts with 80% confidence at supersonic cruise. 
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Sub-scale Modeling issues 


e Cannot directly simulate flight conditions 
— Flight Re = 209 million; boundary layer transition near LE 
— Wind Tunnel Re = 6.4 million; transition is delayed 


e Mixed laminar/turbulent skin friction is too low 
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Sub-scale Modeling Issues: 


Part of the difficulty with wind tunnel testing is that we cannot directly 
simulate flight conditions. 


For example, at flight, the Reynolds number for the TCA configuration is 
209 million. At this Reynolds number, the boundary layer transition from 
laminar to turbulent flow is near the leading edge over most of the wing. On 
the other hand, a much lower Reynolds number is attained in the wind tunnel. 
The Reynolds number for the 1.675% TCA Model in the NASA Langley 
UPWT is equal to 6.4 million. For this low Reynolds number condition, the 
boundary layer transition is delayed over most of the wing. This is 
demonstrated in the flow visualization image below, which shows the extent of 
laminar flow on the model. 


As a result of the mixed laminar/turbulent flow, the skin friction 1s lower 
than that for a flow which is fully-turbulent from the leading edge. 
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Forcing Early Transition 


¢ Standard approach is to add roughness near 
wing leading edge 
¢ This does not completely solve the problem 


~ Transition may not move all the way forward, 
especially in supersonic flow 


— Roughness can create additional drag 


¢ Analysis of wind tunnel drag data requires 
accounting for laminar run & trip drag 
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Forcing Early Transition: 


The standard experimental approach to solve this problem is to add 


roughness (e.g., grit, glass beads, epoxy dots) near the leading edge of the 
wing in order to induce early transition. 


However, this does not always resolve the difficulty. First of all, the 
transition line may not move all the way forward to the trip location, especially 


in supersonic flow. Second, applying roughness to the model can create 
additional drag. 


Therefore, it is important that the analysis of the wind tunnel data takes 
into account all of the laminar run which exists behind the leading edge and 
any additional trip drag that occurs. 
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Not a New Problem 
¢ Braslow (Langley, 1960s) 
— Advocated extrapolating to zero trip height 


¢ Daugherty (Ames, 1970s) 


— Observed constant drag for a range of trip sizes 


¢ Vaucheret (France, 1960s) 


— Applied laminar run correction 


— Corrected drag was constant for small trips 
In all cases, prompt transition was obtained 
for sufficiently large roughness 
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Not a New Problem: 


The whole issue involved with tripping the flow and correcting the wind 
tunnel data is not a new problem in wind tunnel testing. 


Braslow worked at NASA Langley in the 1960s to develop a method to 
estimate trip drag. He advocated using an extrapolation technique. This 
method extrapolates a curve through the drag data as a function of trip height in 
order to estimate the drag for zero trip height. This technique will be further 
examined in Slide 11. 


Daugherty conducted a series of wind tunnel tests at NASA Ames in the 
1970s on the XB-70 model. He observed constant drag over a range of trip 
sizes. Some of his results are presented on the next slide. 


Vaucheret conducted tests on a delta wing/body model in France during the 
1960s. Based on sublimation images, he observed some laminar flow on the 
model. Therefore, he first applied a laminar run correction to his drag data. 
Once he made the correction, he found that the drag was constant for small 
trips. Some data from his work are shown in Slide 10. 


In all cases, it is believed that all of the researchers used airfoil sections with 
sharp leading edges which produced gentle pressure gradients. Therefore, they 
were able to obtain prompt transition for sufficiently large roughness sizes. 
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—- .———___”=. = SeSiaistapaie ee a - — Oras ar 


Jim Daugherty—XB-70 Test 


eS eee mes PReaNenae 
O30 O35 
k, om 


Daugherty, et al. (Ames) 
* M=2.53, Re = 4(10°)/ft, sharp LE sections 
¢ Transition observed at or near the grit strip 


¢ Drag “plateau” observed; zero trip drag for k up to 
about 0.013 in 
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Jim Daugherty—XB-70 Test: 


For further information, please refer to the following references: 


James C. Daugherty, ““Wind-Tunnel/Flight Correlation Study of Aerodynamic 
Characteristics of a Large Flexible Supersonic Cruise Airplane (XB-70-1), 
I—-Wind-Tunnel Tests of a 0.03-Scale Model at Mach Numbers from 0.6 to 
2.53,” NASA Technical Paper 1514, November 1979. 


Henry H. Arnaiz, John B. Peterson, Jr., and James C. Daugherty, “‘ Wind- 
Tunnel/Flight Correlation Study of Aerodynamic Characteristics of a Large 
Flexible Supersonic Cruise Airplane (XB-70-1), I—A Comparison Between 
Characteristics Predicted from Wind-Tunnel Measurements and Those 
Measured in Flight,” NASA Technical Paper 1516, 1979. 


1171 


HSR Annual Airframe Technical Review 


Xavier Vouchatel 


LA KECHERCIE ABKOSESTIARE Ne ise . spuomiine Dateien 1947 
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¢ M= 1.65, Re = 9.8(10°)/m, a = 0°, symmetrical sections 
¢ Trip consisted of a single row of glass beads 
¢ Drag “plateau” found after correcting for laminar flow 
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Xavier Vaucheret: 


For more details of Vaucheret’s work and results for Mach 2.34, please refer to 
the following reference: 


Xavier Vauchert, “Artificial Initiation of the Supersonic Transition,” 
Translation of “Declenchement Artificie] de la Transition en Supersonique,”’ 
La Recherche Aerospatiale, No. 120, pp. 25-32, 1967, NASA TT F-11, 558. 
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Drag vs. Trip Height 


= Analysis per Braslow: 
¢ Plot drag vs. trip height 
¢ Various methods used 


to extrapolate to zero 


= Soe ; ; 
GR TO a BY, eS SH trip height 
Sic SS a > ca TO OR = a =| 5 Gs 

ms ¢ Trip drag is difference 


= ii between measured & 
0.0127 mea tate extrapolated values 
ite eo 
a reset) 6*: Unacceptably large 
shies oe 
0.600 0.005 0.010 0020 range of corrections 


Trip Height (inches) obtained 


How do we resolve this uncertainty? 
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Drag vs Trip Height: 


The Configuration Aerodynamics Test Team attempted to follow 
Braslow’s approach in their initial analysis of the wind tunnel data. The first 
step is to plot the drag as a function of trip height for those values where 
prompt transition is believed to occur. This is shown in the plot on the left. 
Then, a curve is extrapolated through the drag data to zero trip height. The 
trip drag is computed as the difference between the measured drag and the 
extrapolated drag. 


For this data, it is not clear which extrapolation technique is the correct one 
to use. If a quadratic extrapolation is made, the trip drag is equal to 3.3 counts. 
On the other hand, if a linear extrapolation is applied, then the trip drag 
becomes 5.8 counts. The difference between these two methods—2.5 
counts—is unacceptably large. 


At this point, it was important to determine a way to resolve this 
uncertainty in the trip drag. 
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Test 7 


Blunt ft: 


Sharp LE. e 
* Test conducted at Langley UPWT in Jan. 1997 


¢ Collaborators included Mina Cappuccio (Ames) 
and Kevin Mejia (BCAG) 
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Test #1679: 


In order to resolve the uncertainty in estimating trip drag, a more careful 
examination of the results from Test #1679 was undertaken. This test was the 
first entry of the 1.675% Baseline TCA Model 2b in the Unitary Pressure 
Wind Tunnel Test Section #2 at NASA Langley Research Center. The other 
collaborators involved in this wind tunnel test were Mina Cappucio and 
Kevin Mejia. 

It is important to note that the inboard section of this model has a blunt 
leading edge, whereas the outboard section has a sharp leading edge. The 
pressure distribution around the blunt leading edge is different from that 
around the sharp leading edge. The pressure distribution in the outboard 
section with the sharp leading edge is more similar to that obtained in the 
previous historical tests. 
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W/B Trip Drag Study 


¢ Trip height on wing varied 
from k = 0.006” to 0.017” 


Trips were epoxy disks o 
with diameter 0.050”, applied 0.6” from LE 


Fach trip configuration was carefully applied 
and measured for uniformity 


Measure forces & moments at Ma = 2.4 


Obtain series of sublimation images at cruise 
(o% = 3.5°) until stable pattern emerges 
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W/B Trip Drag Study: 


In Test #1679, a trip drag study was conducted on the wing/body 
configuration. The trip height was varied from 0.006 inches to 0.017 inches. 
The trips were epoxy disks with a constant diameter equal to 0.050 inches. The 
disks were applied at 0.6 inches from the wing leading edge. The picture in 
the upper right corner shows an example of the trip disks on the model near the 
leading-edge break. 


Each trip configuration was carefully applied and measured in order to 
maintain uniformity. For each trip configuration, forces and moments were 
measured at the supersonic cruise Mach number, which is equal to 2.4. In 
addition, a series of sublimation images at the cruise angle-of-attack were 
obtained until a stable pattern emerged in the sublimation material. 
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Sublimation Process 


Model is positioned with wings vertical 


Cameras are mounted on each side of test 
section 


Upper & lower surfaces of wing are sprayed 
with sublimation material 


Pictures are taken at 30 sec intervals 


Allows a complete record of the development 
of the transition front 
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Sublimation Process: 


To obtain the sublimation images, the model 1s rotated 90° so that the 
wings are vertical. 


Cameras are mounted on each side of the test section to capture the upper 
and lower surfaces of one wing at the same time. 


The upper and lower surfaces of the wing are sprayed with the sublimation 
material, which is a mixture of fluorene in Genesolv. 


The test section is closed and the tunnel is brought up to the supersonic 
cruise condition. 


Then, pictures are take every 30 seconds until a stable pattern emerges and 
most of the material has sublimed. 


This allows a complete record of the development of the transition 
location. 


1176 


PS IRE TIBOR ORO Oboe tae Rep Ore ve NEE ROD eo ae eee 
, <oppetag 


11 minutes 13 minutes 


AMGoodsell/RA Kennelly Slide 15 





Sublimation Reveals Transition: 


This slide gives examples of sublimation images for the free transition 
configuration at four different times. 


The transition location is apparent in these images because the material 
sublimes more rapidly in areas of turbulent flow. 


After 5 minutes, the initial development of the transition location appears. 
As more time elapses, the transition location becomes more evident. The last 
two images demonstrate that the transition location remains stable with time. 
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Locating Transition 


Transition location was measured as a percentage 
of local chord from the digitized images 
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Locating Transition: 

The sublimation images are digitized and read into a software program like 
Adobe Illustrator to help locate the transition line on the model for each trip 
configuration. 


First, the transition location is drawn on the wing, shown as a thick purple 
line in the picture above. Next, the green guidelines are used to measure the 
transition location as a percentage of the local chord for various spanwise 
Stations. 


This procedure is followed for all of the different trip height 
configurations. 
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w= Free transition, 3.75 counts wwmme ().01 22" dots, 1.19 counts 
m= ().0060" dots, 2.05 counts eee J.0140" dots, 1.12 counts 

0.0081” dots, 1.50 counts <== 0.0170" dots, 0.68 counts 
m=" 0.0102" dots. 1.29 counts 


Outboard: abrupt forward movement 


Inboard: delayed movement 


0.0 
15.8 20.0 35, 38. 35, 30.0 
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Transition Location—Upper: 

This viewgraph gives a summary of the transition location for all of the dot 
heights tested on the upper surface. The thick black line represents the 
location of the free transition line. Note that transition is significantly delayed 
over most of the wing without any trip disks applied to the wing. 


In the outboard section, which has sharp leading-edge airfoil sections, 
there is an abrupt forward movement of the transition location with any 
application of the epoxy disks, even for the smallest disk height. For disk 
heights of 0.010 inches and greater, the transition moves forward to the disk 
location. 


On the other hand, the inboard section behaves differently. The blunt 
leading edge makes it difficult to trip the flow and transition is delayed for all 
disk heights. In fact, for the smallest disk height, the transition location is not 
changed from that for free transition, except in the region near the leading- 
edge break in the planform. Prompt transition is not attained for any of the trip 
configurations, even with the largest trip height of 0.017 inches. 


The dissimilar behavior of the transition line in the inboard and outboard 
sections may be due to differences in the pressure gradient and boundary layer 
development between the two sections. 
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Transition Location—Lower 


=——= Free transition. 2.01 counts comm 0] 22" dots, 0.53 counts 
=== (),0060" dots, 0.75 counts maeese  (.0140" dots, 0.52 counts 

0.0081" dots, 0.64 counts -=—= 0.0170" dots, 0.51 counts 
—= 0.0102" dots, 0.57 counts 


Transition moves promptly 
to trip location 


0.0 
15.8 20.0 : a0, 33. ; Ms 36.9 
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Transition Location—Lower: 


The summary of the transition locations for the different trip heights on the 
lower surface of the wing is shown here. 


The outboard section of the lower surface behaves similarly to that of the 
upper surface. 


However, the transition location in the inboard section displays different 
behavior than on the upper surface. First of all, the free transition location is 
closer to the leading edge than observed on the upper surface. Furthermore, 
the application of most disk heights moves the transition forward to the trip 
location. 
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Methodology 


¢ Laminar run correction 
— Compressible flat-plate skin friction code developed 


~ Compute C, for both mixed laminar/turbulent and 
fully-turbulent flow from LE 


— Difference between these is the laminar run 
correction, AC)),,, 


Once corrected, the data provided guidance 
for extrapolating to zero trip height, 
resolving the uncertainty in AC) 


trip 
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Methodology: 

Once it became evident that a significant amount of laminar flow existed, 
particularly on the inboard, upper, wing surface, it was clear that the wind 
tunnel data needed to be corrected for laminar run before the data could be 
analyzed for trip drag. 


In order to make this laminar run correction, co-author Robert Kennelly 
wrote a FORTRAN code that estimates the skin friction for a compressible 
flow over a flat plate. This code was then used to compute the skin friction for 
a mixed laminar/turbulent flow using the measured transition locations and for 
a flow that is fully-turbulent from the leading edge. The difference between 
these two results is the laminar run correction. 


Once the drag data was corrected and plotted against trip height, the proper 
extrapolation of the data to zero trip height became apparent, which resolves 
the uncertainty in the trip drag. 
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Clutter (Douglas), c. 1959 


¢ Flat plate, compressible, several laminar and turbulent C, 
methods 


¢ Iteratively matches laminar & turbulent momentum 
thicknesses at a specified x-location 


¢ Drag areas (D/q) are summed over all elements; normalized 
by Sref to yield C, for configuration 
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CDf: New Code, Old Method: 





This FORTRAN code has been run on a number of platforms and is available 
through Robert Kennelly at rakennelly @mail.arc.nasa.gov. 
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Drag vs. Trip 
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Drag vs. Tnp Height: 


The reader will now be led through a step-by-step procedure for correcting 
the wind-tunnel data. 


First, the drag data for a=3.5° is obtained by interpolating the 7th-order 
polynomial curves which represent the wind tunnel data. The drag data are 
plotted as a function of trip height for all values tested. The uncorrected drag 
data are shown as blue squares. 


Next, the laminar run correction 1s added to the drag data for each disk 
height. The laminar run correction is largest for free transition and is equal to 
5.8 counts. For a trip height of 0.012 inches, the laminar run correction is 1.7 
counts. The corrected data are represented by the black circles. 


The data, corrected for laminar run, fall into two regions that will be 
described further in the next slide. 
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Drag vs. Trip Height 


¢ Tentatively identify 
“drag plateau” formed 
by first two points 
— Some flattening of the 
curve expected 
— Consistent with 
Braslow’s subsonic data 
& others’ supersonic data 
Linear curve fits 
remaining data well 


— Drag proportional to disk 
frontal area (d * h) 


0.005 0.016 0.015 


Trp Height (in.) 





AMGoodsell/RAKennelly Slide 22 


Drag vs. Trip Height: 


First, a “drag plateau” is defined using data from the first two points—tfree 
transition and 0.006 inches. The plateau is only tentatively identified since 
drag data in between the two endpoints are not available from Test #1679. 
However, some flattening of the drag curve is expected as zero trip height 1s 
approached. In addition, the existence of a drag plateau is consistent with 
Braslow’s analysis of earlier subsonic data and other researchers’ supersonic 
data, including those of Daugherty and Vaucheret. 


The remaining data varies linearly with trip height. For trip heights that 
are approximately equal to the boundary layer height or greater, the drag is 
expected to be proportional to the disk frontal area. Since the diameter 
remains constant in this study and only the height is varying, a linear variation 
of the drag with trip height is appropriate. 
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UAH20 


¢ Difference between 
| : ; linear fit and plateau is 
6.0138 3 : ‘ | trip drag AC rip 


0.0719 


* Trip drag is subtracted 
from data points 


¢ Corrected drags are 
now consistent with 
plateau value 


¢ For this configuration, 
plateau value is the 
fully-corrected cruise 
drag sought 


i005 U.O1U 0.015 
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Final Drag vs. Trip Height: 


The trip drag is computed by taking the difference between the linear fit of 
the data and the value of the drag plateau. 


The trip drag is subtracted from each of the data points that lie above the 
plateau value. The trip drag for the disk height equal to 0.012 inches is 1.9 
counts, which 1s close to the value of the laminar run correction for the same 
disk height. 


The final corrected drag for all disk heights are represented by the solid 
black circles. All of the corrected drag data are now consistent with the value 
of the drag plateau. 


Therefore, for this configuration at the cruise angle-of-attack, the plateau 
value is the fully-corrected cruise drag that the Configuration Aerodynamics 
Test Team has sought. 
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Comparison with CFD Results 
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Comparison with CFD Results: 


In this slide, the final corrected wind tunnel data are compared to three 


Corrected WT data are 
in excellent agreement 
with UPS & CFL3D 
(Baldwin-Lomax) 
OVERFLOW results 
may be low because of 
anomalous Spalart- 
Allmaras behavior near 
leading edge 
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different CFD solutions near the cruise angle-of-attack. 


The corrected wind tunnel data are in excellent agreement with the lift and 
drag results from UPS (Scott Lawrence, NASA Ames Research Center) and 
CFL3D (Chih Fang Shieh, Boeing-Phantom Works, Long Beach). Both of 


these codes used the Baldwin-Lomax turbulence model. 


On the other hand, the drag computed by OVERFLOW (Steve Chaney, 
Boeing Commercial Airplane Group, Seattle) is about two counts lower than 
the wind tunnel data. This is most likely due to the anomalous behavior of the 
Spalart-Allmaras turbulence model which exhibits a “laminar-flow effect” 


near the leading edge of the wing. 


1186 





HSR Annual Airframe Technical Review 


Nacelle Trip Drag Study 





¢ Constant wing trip height: 0.012” 






Three nacelle trip heights (inside & out) 
tested: 0.008”, 0.012”, and 0.014” 


¢ Note: large number of trip disks on nacelles 
(49% of wing total) 


Measure W/B/N/D forces & moments 


Sublimation images obtained at a = 3.5° 
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Nacelle Trip Drag Study: 

In addition to the wing/body trip drag study, a reduced trip drag study was 
conducted on the nacelles. For this nacelle study, the trip disk height on the 
wing was kept constant at 0.012 inches. 


Three different trip heights were applied on the nacelles, both internally 
and externally, at 0.6 inches downstream from the leading edges of the 
nacelles. This results in a large number of disks on the nacelles—the total 
number of trip disks on the nacelles is about 49% of the total number on the 
wing. 

Once again, the forces and moments were measured at the supersonic 
cruise Mach number, 2.4, and a series of sublimation images were obtained at 
the cruise angle-of-attack. 
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Poor visibility of transition location around nacelle 
perimeter requires estimate of laminar run correction 
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Nacelle Sublimation: 


This is an example of a sublimation image on the nacelles for a trip height 
of 0.014 inches. The nacelle transition location is not as easily discernible as 
the wing transition is in these images. In addition, only one of the four 
external surfaces of each nacelle is visible. Therefore, it is difficult to locate 
the transition line on the complete external surfaces of the nacelles due to poor 
visibility in the sublimation region and around the nacelles. As a result, the 
laminar run correction will assume that transition occurs at the trip location. 
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WBND Drag Corrections 


e Laminar run correction 

-Wing: AC),,,,, taken from WB 
study (k = 0.012”) 

—Nacelles: assume transition at 
trip location (four exteriors) 
ACpiamn = 0.000028 

—WT data reduction accounts 
for nacelle interiors 


¢ Trip drag correction 
—Wing: ACh. ipw from WB study 
: : 2 —Nacelle: ACp.i5, scaled from 
3,005 U0) 0.035 0.0265 corresponding WB value 


Nacelle Trip Height (inches) ¢ Three final values agree 
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WBND Drag Corrections: 


The plot on the left shows the drag coefficient for the W/B/N/D 
configuration at the cruise angle-of-attack as a function of the nacelle trip 
height. In addition, the complete laminar run and trip drag corrections for the 
W/B/N/D configuration are shown on this plot. Both the laminar run and trip 
drag corrections consist of a separate wing and a separate nacelle contribution. 


The wing contribution to the laminar run correction is constant for all 
nacelle trip heights and is obtained from the wing/body study for a disk height 
equal to 0.012 inches. The laminar run contribution from the nacelles is based 
on the assumption that transition occurs at the trip location for all trip heights. 
Therefore, the nacelle laminar run contribution is also constant. This laminar 
run correction only takes into account the external surfaces of the nacelles 
because the internal skin friction, which includes the assumption of laminar 
flow between the leading edge and the trip location, is already removed from 
the wind tunnel data. The total laminar run correction (2.0 counts) is added to 
the drag data for each trip height. 


Similar to the laminar run correction, the wing contribution to the trip drag 
is based on the value of the wing/body trip drag for a disk height of 0.012 
inches. The nacelle contribution is computed by scaling the corresponding trip 
drag for the wing/body by 0.49, which accounts for both the internal and 
external dots. The combined trip drag value is subtracted from the drag for 
each disk height. 


Once these corrections are applied, the final values agree with each other. 
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Comparison with CFD Results 


ae =Correctanl WT data (Anes) 
Overtiow/Spatan-Aumarav93 (BCAG) 
A = CELA Batdvtn-Lormax (BPW) 
- 2 ° Once again, corrected WT 
data agree well with 
CFL3D/B-L but differ from 
OVERFLOW/S-A by five 


counts 


¢ The OVERFLOW/S-A 
results have not been 
examined for “laminar flow 
effect” at leading edges of 
nacelles 
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Comparison with CFD results: 

The plot on the left shows the drag polar including both final corrected 
wind tunnel and CFD data for the W/B/N/D configuration near the cruise 
condition. The comparison of the wind tunnel data for the W/B/N/D 
configuration to the lift and drag coefficients computed by CFL3D using the 
Baldwin-Lomax turbulence model show good agreement. 


The OVERFLOW results with the Spalart-Allmaras turbulence model are 
about five drag counts lower than the corrected wind tunnel data. It is believed 
that this under prediction of the drag is, once again, due to the turbulence 
model, but the CFD solutions for the complete configuration have not been 
closely examined. 
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Work in Progress— 
A Closer Look at the B’ Layer 


¢ Why is the inboard, upper surface so 
difficult to trip? 


¢ Why do we see prompt transition on the 
outboard wing? 


¢ Approach: Study boundary layers on flat 
plate and wing (CFD results provided by 
Scott Lawrence/Ames) 
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Work in Progress—A Closer Look at the Boundary Layer: 


The authors are now going to switch directions from focusing on 
methodology to attempting to understand why a modified approach to trip drag 
is required. The main question that comes to mind is “why is the inboard, 
upper surface so difficult to trip?”, whereas, prompt transition on the outboard 
wing is easy to obtain. Our initial approach to answering this question is to 
study the boundary layers on both flat plates and the baseline TCA wing using 
CFD results provided by Scott Lawrence at NASA Ames Research Center. 
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Laminar Boundary Layer Profile 
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Laminar Boundary Layer Profiles: 


The plot on the left gives laminar boundary layer profiles computed by 
UPS at 0.6 inches downstream of the leading edge on a flat plate and on the 
baseline TCA wing at both an inboard and outboard section. The UPS 
solutions for the TCA were obtained at 3.6° angle-of-attack. Also shown on 
the plot, as a solid black line, is the incompressible Blasius solution for a flat 
plate. The heights of the trip disks applied to the wing in Test #1679 are 
indicated on the plot to show the disk heights relative to the boundary layer 
thickness. 


The first observation in comparing the Blasius solution to the UPS solution 
for the flat plate is that compressibility thickens the boundary layer. This 
effect would require a larger disk height to trip the flow. Next, examining the 
differences in the boundary layer profile between the flat plate and TCA wing 
suggests that the gentle, favorable pressure gradient in the outboard section 
tends to slightly thin the boundary layer. This can be seen by comparing the 
solid, blue line with the dashed, green line. The opposite effect is observed for 
the inboard section where the boundary layer 1s thicker than on the flat plate 
(compare the dashed, red line to the solid, blue line). This may be a result of 
the attachment line boundary layer on the blunt leading edge. This increase in 
the boundary layer thickness implies that a larger trip height is required in the 
inboard section, but doesn’t completely explain the why even the largest disk 
height does not cause prompt transition. 
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Mach & Rey, in Boundary Liver 
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Mach and Re, in Boundary Layer: 


The plot on the left gives the distribution of Mach number through the 
laminar boundary layers computed by UPS on a flat plate and the TCA wing. 


This plot shows that the flow at the top of all of the trip disks that were tested 
is supersonic. 


The plot on the right provides the distribution of Reynolds number based 
on local flow properties with the length scale equal to the height off the 
surface. This Reynolds number is being examined because the value of the 
Reynolds number at the trip height, referred to Rex, has been used by 
researchers in the past to predict the required trip height that will cause 
boundary layer transition. This plot shows that the value of the inboard trip 
height Reynolds number is lower than that of the flat plate and the outboard 
section. This further implies that a larger disk height is needed to trip the flow. 
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Summary 
¢ Improved laminar run & trip drag corrections for 
supersonic cruise performance testing 


— Careful flow visualization revealed delayed 
transition, even for the largest disks 


— New code developed for laminar run corrections 
— Corrected data could be analyzed for trip drag 
¢ Despite initial appearances, our data are 
consistent with (most) previous results 


¢ Our approach extends prior work, taking into 
account the effects of blunt leading edges 
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Summary: 

In summary, an improved laminar run and trip drag correction 
methodology for supersonic cruise performance testing was derived. This 
method required more careful analysis of the flow visualization images which 
revealed delayed transition particularly on the inboard upper surface, even for 
the largest trip disks. In addition, a new code was developed to estimate the 
laminar run correction. Once the data were corrected for laminar run, the 
correct approach to the analysis of the trip drag become evident. 


Although the data originally appeared confusing, the corrected data are 
consistent with previous results. Furthermore, the modified approach, which 
was described in this presentation, extends prior historical work by taking into 
account the delayed transition caused by the blunt leading edges. 
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Follow-up “Grit Drag” Test in 4/99 


¢ Repeat portions of Test #1679 


¢ Investigate shorter trip disks to verify that a 
drag plateau exists 


e Study Re & angle-of-attack effects 
¢ Can we improve the present method? 










— Are there more efficient and/or effective trips? 


¢ Intentionally contaminate attachment line 







¢ Different trip configurations on outboard wing 






— Is it really necessary to trip the flow at all? 
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Follow-up Grit Drag Test in 4/99: 

A follow-up grit drag test is planned for April 1999 at NASA Langley’s 
UPWT Test Section #2. First, portions of Test #1679 will be repeated to 
confirm the general trends in the wing/body data which were observed during 
that test and reported in this presentation. Shorter trip disks will be tested on 
the wing and the resulting data will be examined to verify the existence of a 
drag plateau. In addition, some testing will be done to study the effects of 
varying Reynolds number and angle-of-attack on the wing transition location. 


However, the main purpose of the follow-up test will be to try to improve 
the present method, which is time-consuming in both the testing and analysis 
aspects. First, more efficient and/or effect trip mechanisms will be examined. 
For the inboard section, intentionally contaminating the attachment line will be 
attempted. In the outboard region, different trip configurations will be tested. 
Second, it may be possible that for the baseline TCA configuration, it is not 
really necessary to trip the flow at all. If the free transition data corrected for 
laminar run provides the final answer, then testing without the application of 
trip disks is possible as long as other aerodynamic characteristics are not 
adversely affected. 
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This test was done for Configuration Aerodynamics in support of WBS 4.3.1.1, Analysis and 
Testing Methods. 
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There were two objectives for this test. First, was to assess the reasons why there is 
approximately 1.5 drag counts (cts) discrepancy between measured and computed drag 
improvement of the Non-linear Cruise Validation (NCV) over the Technology Concept 
Airplane (TCA) wing body (WB) configurations. The Navier-Stokes (N-S) pre-test 
predictions from Boeing Commercial Airplane Group (BCAG) show 4.5 drag cts of 
improvement for NCV over TCA at a lift coefficient (CL) of 0.1 at Mach 2.4. The pre-test 
predictions from Boeing Phantom Works - Long Beach, BPW-LB, show 3.75 drag cts of 
improvement. BCAG used OVERFLOW and BPW-LB used CFL3D. The first test entry to 
validate the improvement was held at the NASA Langley Research Center (LaRC) UPWT, 
test number 1687. The experimental results showed that the drag improvement was only 2.6 
cts, not accounting for laminar run and trip drag. This is approximately 1.5 cts less than 
predicted computationally. In addition to the low Reynolds Number (RN) test, there was a 
high RN test in the Boeing Supersonic Wind Tunnel (BSWT) of NCV and TCA. BSWT test 
647 showed that the drag improvement of NCV over TCA was also 2.6 cts, but this did 
account for laminar run and trip drag. Every effort needed to be done to assess if the 
improvement measured in LaRC UPWT and BSWT was correct. 


The second objective, once the first objective was met, was to assess the performance 
increment of NCV over TCA accounting for the associated laminar run and trip drag 
corrections in LaRC UPWT. We know that the configurations tested have laminar flow on 
portions of the wing and have trip drag due to the mechanisms used to force the flow to go 
from laminar to turbulent aft of the transition location. 
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This is a photograph of the TCA WB configuration tested. The tripping mechanisms used on 
the wing upper and lower surfaces are 0.05-inch diameter epoxy discs spaced 0.2-inches 
apart, center to center, and are 0.012 inches high. They are placed 0.6 inches aft of the 
leading edge in the streamwise direction on all wing surfaces. The same tripping 
mechanisms, of the same diameter and height, are also placed 1 inch aft of the nose, but 
spaced 0.1 inches apart. 


The TCA WB shown here is the solid wing model known as TCA 2b. The TCA 2b model 
has been tested in the LaRC UPWT three (3) times in test section 2. It also was tested in 
BSWT as part of the high RN test of the NCV model and as part of an Internal Research and 
Development (IR&D) test. 
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TCA 2B WB 
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This is a photograph of the NCV WB configuration tested. The exact same trip configuration 
used on TCA 2b was used for the NCV model. This trip configuration is called the baseline 
trip configuration. Any other trip configuration is considered an alternate trip configuration. 
Alternate trip configurations were also applied to the NCV model to assess if it might be 
more effective to trip the flow and cause the drag improvement to increase. Alternate trip 
configurations were also used to assess the trip drag correction. A similar study of alternate 
trip was done very extensively on TCA 2b during test 1679. 
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This is the drag polar from the BCAG OVERFLOW calculations showing 4.5 cts drag 
improvement of NCV over TCA WB at cruise CL of 0.1 at wind tunnel RN of 4 million/ft at 
Mach 2.4. The Spalart-Allmaras turbulence model was used in OVERFLOW for the N-S 
calculations. 
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This is the drag polar of the BPW-LB CFL3D calculations showing 3.75 cts drag 
improvement of NCV over TCA WB as cruise CL of 0.1 at wind tunnel RN of 4million/ft. 
The Baldwin-Lomax turbulence model was used with CFL3D for the N-S calculations. 
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This figure shows drag versus lift for NCV and TCA WB as measured during test 1687. The 
plot on the left is a close-up of the drag polars around minimum drag for the two 
configurations. The plot on the right is a close-up around cruise. The solid curve is for TCA 
and the dashed curve is for NCV. The curves represent the mean, using a 7” order 
polynomial fit, through all runs within test 1687 of the same configuration. Test 1687 shows 
the drag improvement for NCV over TCA WB is only 2.6 cts. This is approximately 1.5 cts 
less than the calculated pre-test predictions. 
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The approach taken to meet the objectives of the test is to measure forces and moments on 
both the TCA and NCV; get off body flow characteristics using Laser Vapor Screen (LVS) 
for both configurations; get TCA and NCV surface flow characteristics using UV oil and 
mini tufts; and test NCV with alternate trip configurations. 
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Test procedures were carefully looked at and followed throughout the test to make sure the 
most accurate data possible was obtained to reduce errors that might have contributed to the 
measured drag improvement during test 1687. Angle of attack is one of the biggest 
influences on drag. An observation made during the test just prior to this one was that the 
adapter between the sting and strut needed to be replaced. The hardware set up produced an 
error in alpha that was not repeatable. In addition to installing a new adapter, brass buffers 
were added on the retaining area of the sting. The fit between hardware pieces uses a 
cylindrical fit. The buffers allow for the errors in alpha to be reduced when the system 1s 
under load. 


Another component of angle of attack is the droop angle. This is the angle incurred by the 
weight of the model on the hardware support system. This angle gets measured for both the 
upright and inverted runs. The angle should be equal unless there is slop in the system or if 
the fit of the model on the balance is not perfect. For TCA, the upright is different from the 
inverted droop angle. The values of —-2.92° upright and -2.54° inverted were consistent with 
previous tests. For NCV, the upright and inverted droop angles were the same. The value of 
~2.83° upright and inverted was different from test 1687. It is felt that the difference 
between tests is due to the hardware change. While for TCA, the difference between upright 
and inverted is tied up in the fit of the model on the balance. 


An additional component of angle of attack is the balance to model misalignment angle. The 
model was designed such that the balance is installed in the model at a 2° angle down at the 
front of the balance. Models are rarely built such that this angle is met to within 0.001° 
accuracy. Test 1687 used the theoretical value for the balance to model misalignment angle. 
Test 1703 used the NASA Ames Research Center (ARC) QA measurements of —2.0061° for 
TCA and —2.0086° for NCV. What these new values ultimately effected was the computed 
stream angle correction. After test 1703 while QA’ing TCA and NCV, the balance to model 
misalignment angle was measured. There was a larger difference in the measurement of the 
TCA than for the NCV. The TCA LaRC QA measurement was smaller by 0.0166°. The 
NCV LaRC QA measurement was 0.002° smaller. These new measurements, if applied to 
the data, would not affect the data because the stream angle correction would compensate for 
the difference. 


Another large influence to drag is the actual instrument used for measuring forces and 
moments. Balance UT65A was used for all HSR tests in the LaRC UPWT for TCA and 
NCV. The balance was monitored throughout the test. It showed shifts in normal, side, and 
axial force gages. Most of the shifts were within tolerances required by HSR. There was a 
time these shifts were outside of HSR tolerances and it was found that the UV oil damaged 
the moisture proofing on the gages. After this discovery, the balance was removed for the 
gages to be moisture proofed and for the rest of the test the balance was protected during the 
UV oil runs. | 


Another large correction to drag is the pressure acting on the balance by the cavity in the 
aftbody. The aftbody for the NCV configuration has an unusual shape due to the 
optimization process. The TCA aftbody cross section is circular. It was thought that the 
pressures across the aftbody shape might not be constant. Additional pressures were 
measured just behind the balance and just ahead of the fuselage exit at various radial 
locations. 
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The test was heavily supported on and off site. A large amount of flow visualization data 
was acquired. This required good documentation of the results. Images were acquired using 
Hasselblad and high resolution digital cameras. In addition to still images, standard and 
digital video cameras were used for moving images. The images, plots, data, and notes were 
posted on ADAPT daily. ADAPT is a secure website being used within the HSR 
community. No data was ever passed to others within HSR in a non-secure way. 


The reason for posting all this information daily is so that others working on this test that 
couldn’t be on site could participate during the test. Telecons were held daily between 
experimentalists and computational analysts. The team worked together to analyze the 
results and help make any changes necessary in the test plan. 


When ever an interesting feature was detected that hadn’t been looked at as part of the pre- 
test predictions, a computational solution or analysis was done to help the team assess what 


was observed. 


This test could not have been successful as it was without the teamwork of NASA and 
Boeing personnel. 
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The results from the test are presented as outlined for the rest of the presentation. First, the 
short and long term repeatability of the data will be shown. Second, the mini trip drag study 
that was done resulted in a laminar run and trip drag corrections to be made. Finally flow 
visualization images will be shown for the UV oil and LVS runs. 
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This figure shows the short-term repeatability of the lift and pitching moment for TCA. 
These are mean curves through all the repeat runs for this configuration during the test. 
Short-term repeatability was very good for lift and pitching moment. 
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This figure shows the short-term repeatability of drag for TCA. A 7” order polynomial fit 
was put through the repeat runs to get a mean curve. The left plot shows a close-up around 
minimum drag. The right plot shows a close-up around the cruise point. The minimum to 
maximum range of drag at a CL of 0.1 shows a repeatability of less than +0.5 cts. This is 
within the maximum tolerance required by HSR for supersonic cruise performance out of a 
wind tunnel test. 
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This figure shows the long-term repeatability of the lift and pitching moment for TCA. 
These are mean curves through all the repeat runs for this configuration during each test. 
Long-term repeatability was very good for lift and pitching moment between tests 1687 and 
1703. It was slightly worse between tests 1703 and 1679. Test 1679 was the very first time 
the TCA-2b model was tested. The change in angle of attack between tests 17703 and 
previous tests at the CL of 0.1 is 0.0025° for test 1687 and 0.02° for test 1679. 
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This figure shows the long-term repeatability of drag for TCA. These are mean curves 
through all the repeat runs for this configuration using a 7” order polynomial fit through all 
the data during each test. The left plot shows a close-up around minimum drag. The right 
plot shows a close-up around the cruise point. Long term repeatability was very good for 
drag between tests 1687 and 1703. It was worse between tests 1703 and 1679. The change 
in drag between tests 1703 and previous tests at the CL of 0.1 is 0.1 cts less for test 1687 and 
0.7 cts less for test 1679. 
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This figure shows the short-term repeatability of the lift and pitching moment for NCV. 
These are mean curves through all the repeat runs for this configuration during the test. 
Short-term repeatability was very good for lift and pitching moment. __ 
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This figure shows the short-term repeatability of drag for NCV. A 7" order polynomial fit 

was put through the repeat runs to get a mean curve. The left plot shows a close-up around 
minimum drag. The right plot shows a close-up around the cruise point. The minimum to 

maximum range of drag at a CL of 0.1 shows a repeatability of less than +0.25 cts. This is 
within the maximum tolerance required by HSR for supersonic cruise performance out of a 
wind tunnel test. 
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This figure shows the long-term repeatability of the lift and pitching moment for NCV. 
These are mean curves through all the repeat runs for this configuration during each test. 
Long-term repeatability was very good for lift and pitching moment between tests 1687 and 
1703. The change in angle of attack between tests 1703 and 1687 at the CL of 0.1 is 0.025° 


for test 1687. 
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This figure shows the long-term repeatability of drag for NCV. These are mean curves 
through all the repeat runs for this configuration using a 7" order polynomial fit through all 
the data during each test. The left plot shows a close-up around minimum drag. The right 
plot shows a close-up around the cruise point. Long term repeatability was very good for 
drag between tests 1687 and 1703. The change in drag between tests 1703 and 1687 at the 
CL of 0.1 is 0.4 cts less for test 1687. 
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This figure shows the drag improvement for NCV over TCA for test 1703. NCV is 2.3 cts 
less than TCA. Recall test 1687 showed 2.6 cts improvement for NCV over TCA. 
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In an attempt to check if the tripping mechanisms used were the right ones for NCV, CFD 
pre-test calculations were done to assess the boundary layer thickness and the location of the 
baseline trip configuration with respect to the pressure distribution. The long dashed vertical 
lines are the trip discs located 0.6 inch aft of the leading edge in the streamwise direction. 
The dash-dot vertical lines represent the trip discs to be located at 1 inch aft of the leading 
edge. The top two plots are of the pressure distribution on the upper wing surface at two 
spanwise locations close to side of body. It was felt that the baseline location was not 
tripping the flow effectively because they were ahead of the pressure peak. The bottom two 
plots are of the boundary layer and displacement thickness. At 0.6 inch aft of the leading 
edge, the baseline trip discs were 0.012 inches high. This discs are taller than the boundary 
layer thickness. At 1 inch aft of the leading edge, it was felt that 0.015-inch high discs would 
be a good starting point to evaluate the trip effectiveness since this is the height of the 
boundary layer thickness. It was also evident that the trip configuration only needed to be 
replaced on the inboard upper wing surface. The baseline trip configuration was effective 
everywhere else. 
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In addition to the baseline trip configuration and 0.015-inch high discs located 1 inch aft of 
the leading edge, three other trip configurations were investigated as part of the mini trip drag 
study on the NCV model. The trip discs were totally removed from the upper and lower 
wing surfaces to get the data for free transition. The other trip heights investigated on the 
inboard upper wing surface were 0.012 and 0.016 inch high discs. 
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The 0.016 inch high discs on the inboard upper wing surface were the most effective, but still 
did not force the flow to go turbulent just aft of the discs. There still was a small amount of 
laminar run past the disc location. The drag improvement of NCV over TCA using these 
discs on the NCV model increased the improvement only by 0.17 cts from the baseline trip 
configuration for a total drag improvement of 2.57 cts. This data is still not corrected for the 
laminar run and the trip drag. 
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To be able to make the laminar run correction, a series of images were taken during the 
sublimation flow visualization runs. Sublimation is a material that sublimes on the area of 
the wing that is turbulent. The material that remains represent s the areas of laminar run. 
These images are scanned and read into a program like Adobe Illustrator. The transition line 
is laid out on the image and points representing the line as a percentage of the wing chord are 
outputted. This data is used in a boundary layer code that was developed at NASA Ames 
Research center by Robert Kennelly. The code outputs the skin friction drag based on flat 
plate boundary layer theory on a fully turbulent and mixed laminar-turbulent wing flow. The 
difference between these outputs is the drag due to the laminar run that gets added to the total 
uncorrected drag. 


I will be showing the images of some of the sublimation flow visualization runs on the NCV 
model. A trip drag study was done on TCA during test 1679. The results of the laminar run 
and trip drag analysis from the data acquired during test 1679 will be applied to the TCA data 
from this test. 


This image is the transition line laid out for the free transition case on the upper wing surface 


of the NCV model. The points represent the transition line for the boundary layer code to 
compute the laminar run correction. 
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This image is from the sublimation run on the upper wing surface. The trip configuration 
here is 0.016-inch high discs located 1 inch aft of the leading edge inboard and 0.012 inch 
high discs located 0.6 inch aft of the leading edge outboard. To help determining the location 
of transition, two discs were removed on wither side of a disc to create a lone disc in regions 
of the wing. These lone discs help during a sublimation run to determine where transition 
occurs. A 15° turbulent wedge is created where the flow transitions from laminar to 
turbulent flow. These wedges can be seen readily wherever a lone disc exists. This helps the 
transition line to be laid out on the image on the wing surface. The transition line can 
sometimes be difficult to observe due to not being able to distinguish laminar from turbulent 
flow. The best attempt was made to lay out where turbulent flow starts. 
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This is a compilation chart of the transition line for each of the trip configurations tested for 
the upper wing surface. For each transition line, a corresponding laminar run correction is 
listed in the legend. An interesting observation is that the baseline trip was just as effective 
as the alternate trip location except for a small region of the wing near side of body. The 
difference between the baseline and the alternate trip laminar run is less than 0.2 cts. 
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Sublimation flow visualization runs were also done on the lower wing surface for the 2 trip 
configurations tested. This is an image of the transition line laid out for the free transition 
case. 
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This is an image of the transition line laid out for the baseline trip configuration on the lower 
wing surface. 
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This compilation chart show the difference between free transition and the baseline trip 
configuration. Note the amount of laminar run correction made for the lower wing surface 
for each configuration. | 
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This plot shows the uncorrected and laminar run corrected drag at the lift coefficient the 
sublimation runs were acquired. The procedure on how trip drag is determined is described 
in the paper that proceeded this one by Aga Goodsell and Robert Kennelly. Due to this being 
a mini trip drag study, not enough trip configurations were tested to be able to determine 
where the drag plateau occurs. An estimation was made that the plateau occurs at zero trip 
height and that the rest of the data needed to cross the plateau at a trip height of 0.007 inches. 
These estimations came about from the trip drag study done on TCA during test 1679. A 
linear fit is put through the data for the configurations tested with trip discs. The difference 
between the linear fit and the plateau is the trip drag correction. 
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This plot shows the drag corrections made for laminar run and trip drag for each 
configuration tested. These corrections are applied to the whole drag polar. The laminar run 
correction is added, while the trip drag correction is subtracted. 
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After making these corrections, the drag improvement of NCV over TCA only improved by 
0.2 cts more than that determined using the uncorrected data. This applies to all tests and all 
trip configurations. For the baseline trip configuration the improvement of NCV over TCA 
is 2.8 cts. For test 1703, the improvement is 2.5 cts. The difference between the 2 tests is 
wrapped up in the long-term repeatability. For the alternate trip location and the highest 
height, the improvement is 2.8 cts. The difference between the baseline and the alternate trip 
configuration is wrapped up in the cumulative error in laminar run, trip drag, and data 


repeatability. 
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This is a Summary image of the UV oil run done on TCA at the supersonic cruise angle of 
attack. The upper most image is of the upper wing surface. Just under the upper wing is the 
lower wing image from the UV oil run. At the bottom left is side of body near the wing apex 
and on the right is looking at the forebody upper surface. Hasselblad cameras were mounted 
on both windows. This image makes up 2 runs. One run was with the model mounted wings 
vertically to get the upper surfaces. The second run was with the model mounted wings 
horizontally to get side of body. The images depicts the surface flow on the model. The 
flow is very clean without any separation. 
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This is a compilation of images from UV oil runs on the NCV model at the supersonic cruise 
angle of attack. Here you will notice that the upper wing surface has some distinct features 
representing vortices. There seems to be two traces of the vortices on the wing. One goes 
streamwise over the inboard nacelle location and the other over the outboard nacelle location. 
The one that passes over the inboard nacelle location seems to be emanating from the wing 
body junction. The one that passes over the outboard nacelle location seems to come from 
the same region but moves span wise as it goes downstream. The vortices will be evident 
when the LVS images are seen. 


The flow over the lower wing surface, side of body, and upper forebody look clean. There is 
no indication of separation in any of these regions. 
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For every UV oil run, images for a maximum of three (3) angles of attack could be obtained. 
This figure shows the UV oil results for the cruise and two (2) higher angles of attack for the 
NCV model. The two (2) vortices that seem to pass over the upper wing surface become 
more evident at the higher angles of attack. They become stronger and cause flow separation 
on the wing surface. 
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NCV Wing UV Oil for Different Alphas 
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Alpha=8° 


Alpha=6° 


Alpha=4.5° 


This figure show the UV oil results on the NCV model at the same three (3) angles of attack. 
The flow turns more and the attachment line moves as the angle of attack increases. 
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Post test, a comparison was made between the UV oil images and the CFD calculations that 
were done at NASA Ames Research Center by Scott Lawrence using UPS. The top figure 
compares the fully turbulent solution from UPS with the UV oil image. The bottom figure 
compares it with a fully laminar solution. It was determined from this that the original 
thought that the upper wing surface had laminar flow on it was incorrect. The flow is 
turbulent except ahead of the trip mechanisms and slightly behind. 
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The UPS solution was also compared to the UV oil images for side of body. It was feared 
that the flow over the NCV forebody was laminar. It was hard to see during the sublimation 
runs. The top figure compares the fully laminar UPS solution with the UV oil image. The 
bottom two (2) figures compare fully turbulent UPS solutions with the UV oil image. The 
turbulence models being compared are Spalart-Allmaras and Baldwin-Lomax. The flow over 
the forebody is more like the fully turbulent UPS solutions. 
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This figure shows the set up used for the Laser Vapor Screen runs. The top left image shows 
the Spectra-Physics laser used. The right image shows a light box attached to the tunnel 
window. The model was translated axially so that the laser light sheet was at a different 
location on the model for tracking any shock and vortices. A Hasselblad camera was 
mounted in a housing to the tunnel ceiling and was remotely triggered in the control room to 
take the still image. A mini video camera was in a housing mounted to the hardware holding 
onto the model. The two (2) cameras can be seen in the bottom left image. 
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Laser Vapor Screen (LVS) Setup 


ANA Teale. . om 
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The model was painted black with reflective targets placed 4 inches apart from the aftbody 
forward. These targets were used to stop the axial movement of the model when the laser 
light sheet passed to take the still images. A video was recorded from the mini video camera 
signal of the entire LVS run. Any shocks and vortices were tracked for every angle of attack 
during the LVS run. 
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During the LVS run water was injected in the tunnel to create condensation to help better 
pick out shocks and vortices. This figure compares the LVS image between TCA and NCV 
at airplane station 2671 inches. The very dark area is the wing behind the light sheet. Above 
the light sheet close to the leading edge, a compression shock can be seen sitting on the wing. 
There is also a horseshoe vortex at side of body on the NCV model. It is hard to see it in this 
figure but the video clearly shows it. The TCA does not have a horseshoe vortex as strong as 
on the NCV as observed from the mini video camera. 
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This figure compares the LVS run for TCA and NCV at the aftbody location, which is 
downstream of the wing trailing edge. This figure helps show that the NCV model has these 
two (2) vortices that pass over the upper wing surface. The NCV wake shows traces of two 
(2) vortices. The TCA wake is clean. 
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This figure compares TCA and NCV LVS results at an angle of attack of 1.5° above the 
cruise angle of attack for each configuration at airplane station 2671 inches. For TCA, the 
image is for 5°, and for NCV the image is for 6°. The dark spots above the light sheet for 
NCV are the horseshoe vortex and the vortex that passes over the inboard nacelle location. 
These cannot be seen for TCA. 
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This figure shows the wake behind the wing trailing edge for TCA and NCV at 1.5° higher 
than the cruise angle of attack. The NCV wake clearly show three (3) vortices: a horseshoe 
vortex at side of body, a vortex that passes over the inboard nacelle location, and a vortex 
that passes over the outboard nacelle location. The TCA wake is just starting to show very 
weak vortices over the wing. 
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The forces and moments acquired on the NCV and TCA models are as good as can be with 
the angle of attack and balance instrumentation used during the test. Both the low and high 
RN tests did not measure the calculated drag improvement of NCV over TCA. Both tests 
show that NCV is approximately 2.8 cts better than TCA. 


Both UV oil and LVS are good flow visualization tools for acquiring surface and off body 
flow characteristics. The two techniques go together well. LVS backs up what is seen in the 
UV oil images. 


The only items that are needed to close-out the experimental portion of the NCV testing are 


to get surface pressure measurements, CFD solutions of as-measured (QA) geometry, and 
validate effect of balance bore angle measurements by LaRC on experimental results. 
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Aft Body Closure - Predicted 
Strut Effects at M=2.4 


John Ek. Lamar, NASA LaRC 
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HSR Annual Airframe Technical 
Review 


Anaheim Marriott, Anaheim, CA 
February 8-12, 1999 





This paper reports the predicted M = 2.4 strut-interference effects on a closed 
aft body with empennage for the TCA baseline model. The strut mounting 
technique was needed in order to assess the impact of aft-end shaping, i.e. open 
for a sting or closed to better represent a flight vehicle. However, this technique 
can potentially lead to unanticipated effects that are measured on the aft body. 
Therefore, a set of computations were performed in order to examine the closed 
aft body with and without strut present, at both zero and non-zero angles of 
sideshp (AOS). The work was divided into a computational task performed by 
Javier A. Garrriz, using an inviscid (Euler) solver, and a monitoring/reporting 
task done by John E. Lamar. All this work was performed during FY98 at the 
NASA Langley Research Center. 
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Goals Demonstrate Significant L/Dmax Gains 


Objectives Robust Analysis / Realistic Aerodynamic Efficient Engine / 
Testing Methods Design Optimization Airframe Integration 


ne) en ae EN mecvoeeites 
Challenges [Validation | Viscous Effects } Multipoint Conditions 


Analytic Methods | Design iq Test Programs 
pperecies and Applications Develo pment and Techniques 


Methods Down Select Nacelle /Diverter | [1 wr Database 
Design Integration 


Program [-|_Viscous Drag Prediction 
ee — Tech. Baseline 
Cruise Point Optimization Development High Re No. Testing 
Mull-Point Optimization Aero S&C 
S&C CFD Predictions evae ment 





The work reported in this paper is associated with the highlighted boxes. 
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Outline 


¢ CFD Objectives and Geometrical Modeling 


¢ Supersonic Results 


Conclusions 


Recommendations 





This presentation has four parts. The first deals with the objectives of the CFD 
study and details the geometrical modeling employed. The second part shows 
the resulting pressures and forces/moments obtained at M = 2.4 with the strut- 
off and -on. The last two parts summarize the paper with some conclusions and 
provides recommendations for further studies. 
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CFD Objectives and 
Geometrical Modeling 


¢ Estimate strut interference on TCA airplane 
model at M=2.4, 0.6, 0.9 and 1.2 for AOA 
4and AOS of 0 and -3 using an Euler solver 
(USM3Dns in its inviscid mode) 


¢ Complete airplane - tip to tip - and 
supersonic strut modeled 


¢ Strut added to airplane at AOS of 0 and -3 





The general goal of this task 1s to be able to estimate the strut interference effects on the 
TCA airplane model at both transonic and supersonic speeds. To date only the M = 2.4 
com putations have been completed at an angle-of-attack (AOA) = 4 degrees and at 
AOS = 0 and -3 degrees. These have been made using the USM3Dns code in its Euler 
or inviscid mode. An overall description of the code follows: 


USM3 Dns 1s a tetrahedral-based Euler and Navier-Stokes flow solver. Spatial 
discretization is accomplished by a cell-centered finite-volume formulation using an 
accurate linear reconstruction scheme and upwind flux differencing. Time is advanced 
by an implicit backward-Euler time-stepping scheme. Flow turbulence effects are 
modeled by the Spalart-Allmaras one-e quation model, which is coupled with a wall 
function to reduce the number of cells in the sublayer region of the boundary layer. 


Since the study involved sideslip, it was necessary to model the entire configuration and 
not take advantage of symmetry. Even though this was not required at AOS = 0 
degrees, for consistency full configuration modeling was maintained for all 

com putational situations. 


In addition, the supersonic strut was modeled, including the region where it abutted the 
lower surface of the model. All together about 1.3 million cells were used to model the 
com plete airplane and strut combination. 


During the computations, the model geometry was maintained at AOS = 0 degrees 
regardless of the onset flow. For AOS = -3 degrees the boundary conditions for the 
flow solution were changed on the model and the relative alignment angle (as viewed 
from the top) between the model and the strut was altered by moving the strut only the 3 
degrees. This was done because the strut surface was represented by a fewer number of 
cells than the model. 
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Perspective Top View of TCA 





This figures shows a top view of the complete TCA aircraft model from a 3/4 
front perspective. The different colors on the various parts of the model are 
only to assist in keeping track of the geometrical elements during the grid 
generation portion of this task. 
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Top View of TCA - UPWT 
Strut at AOS of 0 





This figure shows the arrangement of the complete TCA model and the 
supersonic strut from a top view at an AOS = 0 degrees. Again the various 
colors shown on the different surfaces are only used to keep track of the 
geometrical elements during the grid creation process. 
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Side View of TCA - UPWT 
Strut at AOS of 0 





Side view of TCA model mounted atop the supersonic strut at AOS = 0 degrees. 
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Perspective View of TCA - 
UPWT Strut at AOS of 0 





This figure is a lower surface view from a 3/4 front perspective of the TCA 
model mounted atop the supersonic strut at AOS = 0 degrees. Some details of 
the strut system, both in terms of its cross sectional change with longitudinal 
distance and the face of the strut where attaches to the tunnel support system, 
can be seen here. 
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Alt. Perspective view of TCA - 
UPWT Strut at AOS of 0 





This figure is another lower surface view of the TCA model and supersonic 
strut system at AOS = 0 degrees, but this time looking forward from underneath 
the model and with a 3/4 rear perspective. 
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Top View of TCA - UPWT 
Strut at AOS of -3 





This figure shows the arrangement of the complete TCA model and the 
supersonic strut from a top view at an AOS = -3 degrees. Again the various 
colors shown on the different surfaces are only used to keep track of the 
geometrical elements during the grid creation process. 
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Perspective View TCA - UPWT 
Strut Juncture at AOS of -3 





This figure is a lower surface view from a 3/4 front perspective of the juncture 
between the TCA model mounted and the supersonic strut at AOS = -3 degrees. 
Note that at the top of the strut leading-edge there 1s an air gap between the strut 
and the model which does occur because of the 3 degree misalignment. This 
gap was modeled. 
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Supersonic Results 


e Surface Pressure Contours 


¢ Forces/Moments 


e« Cross Sections and Pressures 





The figures that follow present the M = 2.4 aerodynamic results for AOA = 4 
degrees at AOS = 0 and -3 degrees with the strut-off and -on. These figures 
include eight surface pressure contours, four for the upper and four for the 
lower; the overall forces/moments for the entire configuration, as well as the 
forces/moments for the aft body alone; and cross-sectional plots of local 
pressures and associated wing/body geometry at fixed Fuselage Stations. The 
figure presentation follows this order. 
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TCA CONFIGURATION WITHOUT UPWT STRUT 
M=2.4,0=4,f=0° 


USM3D Inviscid Flow Solution 
Upper Surface Cp 
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This figure shows the upper surface pressure coefficient contours for the TCA 
confi guration without the supersonic strut at M = 2.4, AOA = 4 degrees, and 
AOS = 0 degrees. The most noticeable phenomenon is that even at AOS = 0 
degrees the contours are not exactly symmetrical. This 1s due to the 
unstructured CFD code (USM3D) model img the entire configuration and not 
taking advantage of geometrical symmetry. Consequently, one should expect 
small lateral side-force and moments at this attitude. 


Along the inboard leading-edge and trailing aft over the outboard wing, a region 
of higher suction pressure is present. 
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TCA CONFIGURATION WITH UPWT STRUT 
M=24,0a=4,{8=0 


USM3D inviscid Flow Solution 
Upper Surface Cp 





This figure shows the upper surface pressure coefficient contours for the TCA 
confi guration with the supersonic strut at M = 2.4, AOA = 4 degrees, and AOS 
= 0 degrees. In comparison with the previous figure (strut-off), there are no 
striking differences seen on the upper surface associated with the strut, as would 
be expected under these conditions. 
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TCA CONFIGURATION WITHOUT UPWT STRUT 
M=24,a=4,8=0 
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USM3D Inviscid Flow Solution 
Lower Surface Cp 
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This figure shows the lower surface pressure coefficient contours for the TCA 
confi guration without the supersonic strut at M = 2.4, AOA = 4 degrees, and 
AOS = 0 degrees. Just as for the upper surface, these contours are not exactly 


symmetrical even at AOS = 0 degrees. The most interesting features here are 
the shock-on-shock interactions between the nacelles, the relatively large 
suction values of pressure coefficient reached inboard of the inner nacelles, due 
to flow expansion, and the larger positive pressure coefficients near the outer- 
wing leading edge. 
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TCA CONFIGURATION WITH UPWT STRUT 
M=2.4,0=4,8=0° 


USM3D Inviscid Flow Solution 
Lower Surface Cp 





This figure shows the lower surface pressure coefficient contours for the TCA 
confi guration with the supersonic strut at M = 2.4, AOA = 4 degrees, and AOS 
= 0 degrees. In comparison with the previous figure (strut-off), the nacelle 
shocks and outer-wing leading edge pressure coefficients appear the same. 
However, there are differences in the pressure coefficients in the streamwise 
region between the fuselage/wing juncture and the start of the wing outer panel. 
The presence of the strut is to cause this region of the wing (an inverted V- 
shape just inboard-and-downstream of the inboard leading-edges) to have a 
larger positive value than without the strut. Moreover, along the centerline and 
downstream of the wing/fuselage juncture -- where the strut footprint would be 
located (see figure 25) -- there is a region of higher suction pressure coefficient 
indicated on the model. This is unexpected result for a no-flow situation and 
therefore the integration of the pressure underneath the strut footprint is not 
performed. Another reason for omitting these force/moment contributions to 
the reported values is given on later on figure 21. 
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TCA CONFIGURATION WITHOUT UPWT STRUT 
M=2.4,0=4°,B=-% 


USMS5D Inviscid Flow Solution 
Upper Surface Cy 





This figure shows the upper surface pressure coefficient contours for the TCA 
confi guration without the supersonic strut at M = 2.4, AOA = 4 degrees, and 
AOS = -3 degrees. The effect of sideslip is clearly noted on the inboard wing 
where the higher suction region on the windward wing moves inboard and on 
the leeward wing moves outboard. Moreover, asymmetries are noted in local 
pressures over the entire length of the model, including higher suctions on the 
leeward portion of the horizontal tail. 
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TCA CONFIGURATION WITH UPWT STRUT 
M = 2.4,0=4°, 8 =-3° 


USMSD Inviscid Flow Solution 
Upper Surface Cp 





This figure shows the upper surface pressure coefficient contours for the TCA 
confi guration with the supersonic strut at M = 2.4, AOA = 4 degrees, and AOS 
= -3 degrees. In comparison with the previous figure (strut-off), there are no 
striking differences seen on the upper surface associated with the strut, as would 
be expected under these conditions. 
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TCA CONFIGURATION WITHOUT UPWT STRUT 
M=2.4,0=458=-3° 


USM3D inviscid Flow Solution 
Lower Surface Cp 
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This figure shows the lower surface pressure coefficient contours for the TCA 
confi guration without the supersonic strut at M = 2.4, AOA = 4 degrees, and 
AOS = -3 degrees. In comparison with figure 15 at AOS = 0 degree, the 

asy mmetries associated with sideslip are present all along the configuration, but 
especially noticeable are the changes to the nacelle shocks and the pressures on 
the lifting surfaces. 
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TCA CONFIGURATION WITH UPWT STRUT 
M=2.4,0=4°B=-3° 


USM3D inviscid Flow Solution 
Lower Surface Cp 





This figure shows the lower surface pressure coefficient contours for the TCA 
confi guration with the supersonic strut at M = 2.4, AOA = 4 degrees, and AOS 
= -3 degrees. In comparison with figure 19 (strut-off), the nacelle shocks and 
outer-wing leading edge pressure coefficients appear the same. However, there 
are differences in the pressure coefficients in the streamwise region between the 
fuselage/wing juncture and the start of the wing outer panel. In particular, the 
inverted V-shape noted in figure 16 at AOS = 0 degrees is more-or-less retained 
at sideslip. Hence, the strut tends to have a mitigating effect on the pressures in 
its immediate vicinity. 
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TCA Total Force Coefficients 
Estimated by Euler Solver (USM3D) 


M = 2.4, AOA = 4° 


AOS = 0° 
Normal Axial 
Without strut 0.11392 0.00104 
With strut 0.11363 0.00145 


AOS = 3° 


Normal Axial Side 
Without strut 0.11277 0.00097 0.00035 
With strut 0.11316 0.00147 0.01514 





[Note that the forces/moments reported in figures 21 to 24 for the TCA airplane 
model do not contain those generated over the area associated with the strut 
footprint (see figure 25). This is done to keep the comparisons consistent with 
the minimum exposed model-surface-area, 1.e. strut on. ] 


The strut effect on the normal- and side-force coefficients at both AOS of 0 and 
-3 degrees is seen to be small, and its effect on the axial force coefficient at 
either AOS value amounts to a “pressure drag” increase of about ~SxE-04 ( ~5 
drag counts). 

(Note that the side force coefficient should be zero at AOS = 0 degrees and it 
isn’t due to the unstructured grid modeling of the complete configuration.) 
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TCA Total Moment Coefficients 
Estimated by Euler Solver (USM3D) 


M = 2.4, AOA = 4° 


AOS = 0° 
Pitch Roll 
Without strut -0.21010 -0.00229 
With strut 0.00203 -0.00001 


AOS = 3° 


Pitch Roll 
Without strut -0.02056 -0.04442 
With strut 0.00155 0.00157 





At both AOS of 0 and -3 degrees the strut effect on the rolling- and yawing- 
moment coefficients is seen to be small*; however, there is a destabilizing 
effect on the pitching moment coefficient to the extent that it goes from stable 
to slightly unstable. 


* Note they should be zero as should all lateral moments values at AOS = 0 degrees. 
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TCA Total and Aftbody Force Coefficients 
Estimated by Euler Solver (USM3D) 


M = 2.4, AOA = 4°, AOS = 0° 
Total Configuration 


Normal Axial 
Without strut 0.11392 0.00104 
With strut 0.11363 0.00145 


Aftbody Only 


Normal Axial 
Without strut 0.00075 0.00028 
With strut 0.00047 0.00032 





The strut effect at AOS of 0 degrees on the total TCA configuration force 
coefficients are repeated from figure 21 for comparison with those for the aft 
body only. For both it is noted that the changes in aft body force coefficients 
mirror those of the total configuration. In particular, with the strut present there 
is less normal force and more axial force; moreover the increase in the latter or 
“pressure drag” amounts to about one-tenth (~5E-05 or ~0.5 drag counts) of the 
total increase (~5E-04 or ~5 drag counts). 
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TCA Total and Aftbody Moments Coefficients 
Estimated by Euler Solver (USM3D) 


M = 2.4, AOA = 4°, AOS = 0° 


Total Configuration 
Pitch Roll 
Without strut -0.21010 -0.00229 
With strut 0.00203 -0.00001 


Aftbody Only 


Pitch Roll 
Without strut -0.00079 0.0 
With strut -0.00046 0.0 





The strut effect at AOS of 0 degrees on the total moments coefficients are 
repeated from figure 22 for comparison with those for the aft body only. The 
only consistent strut effect noted between the two sets of data is the slightly 
destabilizing effect on pitching moment coefficient. The other aft body moment 
coefficients should all be zero but are not due to the unstructured modeling of 
the complete configuration. 


1496 


TCA Pertinent Locations 


Cp Computation Stations 


To «UPWT Strut °Foorprin’ 
@ Moment Reference Pom (Pianar Location) 





This figure shows the lower surface of the TCA airplane model, the footprint of 
the strut at AOS = 0 degrees, the moment reference point and the thirteen 
computational stations where the local pressure coefficients were computed. 
Not all thirteen locations of pressure data are presented, but enough of them are 
shown in order to provide understanding as to how the reported forces/moments 
are developed. Stations 2 to 7 bound the strut footprint on the model. 
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TCA Configuration 
_ Grid Cross—Section‘and 
Surface Cp Comparisons at Station 1 
 M=2.4,a=4.0° B =0.0° 
USMSD Inviscid Flow Simulation 





[Note for figures 26 to 37 that the view is looking upstream with the right hand 
side having the positive y values. ] 

This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 1, M = 2.4, AOA = 4 degrees and AOS = 
Q degrees. This station was chosen because it is ahead of the strut, hence no 
effect is expected here and thus it can be used for validating the solutions. The 
pressure data show a negligible effect of the strut . 
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TCA Configuration 
Grid Gross—Section and 
Surface Cp Comparisons at Station 3 
M = 2.4,0.=4.0° §:=0.0° — 
USMSD Inviscid Flow Simulation 


0.2 


0.1 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 3, M = 2.4, AOA = 4 degrees and AOS = 
0 degrees. This is the first station on the strut and shows there to be more 
positive pressure on the lower surface with the strut present, hence more normal 
force. Since this station is ahead of the moment reference point, it will 
contribute a more nose-up moment and be destabilizing. 
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TCA Configuration 
Grid Cross-Section and 
Surface Cp Comparisons at Station 4 
M=2.4,0=4.0 §.=0.0 
USM3D inviscid Flow. Simulation 


0.2 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 4, M = 2.4, AOA = 4 degrees and AOS = 
Q degrees. This station also shows there to be more positive pressure on the 
lower surface with the strut present, hence more normal force. Since this station 
is also ahead of the moment reference point, it will contribute a more nose-up 
moment and be destabilizing. 
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TCA Configuration 
Grid Cross-Section and 
Surface Gp. Comparisons at Station 7 
M = 2.4,0.=4.0°,8.=0.0° 
USM3D Inviscid Flow Simulation 


0.2 


-200. 
—400 -200 6 200 400 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 7, M = 2.4, AOA = 4 degrees and AOS = 
0 degrees. This station shows that with the strut present there are two suction 
pressure spikes on the lower surface near y = 0. However, the strut on pressures 
are more positive outboard of y = +/-140, hence the resulting normal force 
could be slightly less or the same-as the no-strut configuration. Since this 
station is still ahead of the moment reference point, the strut effect could be 
slightly stabilizing or produce no change. 
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TGA Configu ration 
Grid Gri section and 
Surface Cp Comparisons at Station 9 
M=24, a=4.0, §.=0:0° 
USM3D Inviscid Flow Simulation 


0.2 


~450 * : 
—600 -300 0 300 600 


woe , | 
-§00 -300 0 300 600 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 9, M = 2.4, AOA = 4 degrees and AOS = 
0 degrees. This station shows that with the strut present there is less positive 
pressure on the lower surface, hence less normal force. Since this station is 
behind the moment reference point, it will contribute a smaller nose-down 
moment and therefore be destabilizing. 


1502 


TCA Configuration 
Gnd Cross-Section and 
Surfaca Cp Comparisons at Station 12 
M= 2.4, @=4.0,, 8.=0.0° 
USM3D Inviscid Flow Simulation 


0.2 


150 : 
—-300 -150 0 ~0.2 
y ~300 ~-150 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 12, M = 2.4, AOA = 4 degrees and AOS 
= 0 degrees. No discernable differences are noted in the pressure coefficients so 
no incremental contribution to the normal force or pitching moment from the 
strut is expected on the empennage. 
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TCA Configuration 
Grid Cross-Section and 
Surface Cp Comparisons at Station 1 
M=2.4,.a=40°8 =-3:0° 
USM3D Inviscid Flow Simulation 


0.2 


Note: Sting should read Strut 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 1, M = 2.4, AOA = 4 degrees and AOS = 
-3 degrees. This station was chosen because it is ahead of the strut, hence no 
effect is expected here and thus it can be used for validating the solutions. The 
pressure data show a negligible effect of the strut. 
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TCA Con figuration 
Grid Cross—Section and 
Surface Cp Comparisons at Station 3 
M=2.4,a=40, B.=-3.0° 
USM3D Inviscid Flow Simu lation 


0.2 


Note: Sting should read Strut 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 3, M = 2.4, AOA = 4 degrees and AOS = 
-3 degrees. This is the first station on the strut and shows there to be more 
positive pressure on the lower surface with the strut present, hence more normal 
force. Since this station is ahead of the moment reference point, it will 
contribute a more nose-up moment and be destabilizing. 


[Note that there will be small side force coefficient and rollmg- and yawing- 
moment coefficient differences associated with the integration of the two 
different lateral, delta-pressure-coefficient distributions over the cross-sectional 


geometry for figures 33 to 37.] 
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TCA ‘Configu ration — 
Grid-C Section and 
Surface. Cp Comparisons at Station 4 
M = 2: 4, 0.=4.0 , B=-3.0° 
nviscid Flow Simulation 


0.2 


-150 -75 O BD 
¥ 


Note: Sung should read Strut 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 4, M = 2.4, AOA = 4 degrees and AOS = 
-3 degrees. This station also shows there to be more positive pressure on the 
lower surface with the strut present, hence more normal force. Since this station 
is also ahead of the moment reference point, it will contribute a more nose-up 
moment and be destabilizing. 
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TCA Configuration | 
Gri —Section and. 
Surface Cp Comparisons at Station 7 
M=2.4,a=4.0 °8.=-3.0° 
USMS3D Inviscid Flow Siralation 


Note: Sting should read Strut 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 7, M = 2.4, AOA = 4 degrees and AOS = 
-3 degrees. This station shows that with the strut present there are two suction 
pressure spikes on the lower surface near y = 0. However, the strut on pressures 
are more positive outboard of y = -100 and +150, hence the resulting normal 
force could be slightly less or the same-as the no-strut configuration. Since this 
station is still ahead of the moment reference point, the strut effect could be 
slightly stabilizing or produce no change. 
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TCA Contiguration 
é Gr Gruss ecipha 6 
urface Cp Comparisons at:Station 
M-= 2:4,0=4.0;8.=-3.0" 
USM3D inviscid Flow Simulation 


0.2 


0.3 , 
—-600 -300 0 300 600 
y 


Note: Sting should read Suut 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 9, M = 2.4, AOA = 4 degrees and AOS = 
-3 degrees. This station shows that with the strut present there 1s less positive 
pressure on the lower surface, hence less normal force. Since this station is 
behind the moment reference point, it will contribute a smaller nose-down 
moment and therefore be destabilizing. 
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TCA Configuration 
Grid Cross-Section and 
Surface Cp Comparisons at Station 12 
= <.4, o.=.4. | 1 = 


USM3D Inviscid Flow Simulation 
0.2 
OL 


0.0 
~C,, 


109 
-300 -150 0 
y 


Note: Sting should read Strut 





This figure shows the grid cross-section and the USM3Dns inviscid solution for 
surface pressure coefficients at station 12, M = 2.4, AOA = 4 degrees and AOS 
= -3 degrees. No discernable differences are noted in the pressure coefficients 
so no incremental contribution to the normal force or pitching moment from the 
strut is expected on the empennage. 
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Conclusions 


Inviscid (Euler) strut effects at M=2.4 on the 
TCA configuration and aft body are: 


- Negligible for normal force coefficient 


- Increase pressure drag (axial force) coefficient 


- Slightly pitch destabilizing 





The USM3Dns code has been used in its inviscid (Euler) mode to provide 
solutions for the TCA model in free air and in the presence of a support strut at 
M = 2.4 at both Angles Of Sideslip (AOS) of 0 and -3 degrees. The strut effects 
on the TCA confi guration and aft body were determined to be negligible on 
normal force, to mcrease the pressure drag by ~5 counts on the TCA 

confi guration and ~0.5 counts on the aft body, and be slightly pitch 
destabilizing on both. 
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Recommendations 


¢ Compare Euler and Navier-Stokes solutions 
over the TCA aft body for particle traces, 
Cp, and force/moment 

¢ Do transonic solutions of transonic strut 


and compete airplane prior to testing in 
16FTT in FY99 





It is recommended that viscous (Navier-Stokes) solutions be generated for 
comparison with those shown herein. The resulting comparisons should include 
particle traces, pressure coefficients, and forces/moments over the aft body. 


It is also recommended that prior to the transonic wind-tunnel testing of the 
TCA model in the 16FTT during Spring 1999, transonic solutions, of at least the 
Euler type, be available to provide understanding of the anticipated strut- 
interference effects. Since the transonic strut is shaped differently than the 
supersonic one, there is a need for some new geometric modeling to be 
completed in order for this task to be preformed. 
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